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NOTES  OH  A THEORY  OF  SPINNING  SHELL 


ABSTRACT 


The  theory  expounded  here  Is  based  esiuly  on  that  of  FOeler, 
Gallop,  Lock  and  Richmond*  except  that  the  coordinate  system  is 
attached  to  the  actual  trajectory  instead  of  the  particle  tra- 
jectory. They  are  followed  in  their  fortunate  choice  of  symbols, 
especially  for  the  damping  factors. 

ihe  equations  for  the  complex  year  are  developed  on  the  usual 
linear  assumption  and  tbsir  solutions  obtained.  Ths  dynssical 
stabjJLity  of  ?b?3I  1*  dlaeuieaii.  It  is  here  pointed  out  that  , 
for  practical  purposes,  it  is  not  essential  that  a shell  be 
dynamically  stable  at  all  points  of  Its  trajectory.  A rc«*&urlng 
comparison  is  made  between  the  results  of  this  theory  and  those 
of  FGLR  and  Kelley  and  McShaee. 


\ 

Phil.  Trane.  Roy.  Soc<=  A,  voi.  221,  pp.  295-5 8?  (1920)* 


>' 


DmIkitiok  of  btmbcls 


I 


r 


a.  Dynamical  constants  and.  variables. 

b.  Geometric  m&gni’-ud^s  and  vectors? 
d.  Aerodynamic  forces. 

d.  Damping  factors. 


a.  Dynamical  constants  and  variables: 

A Ariel  docent  of  inertia  of  ths  shell 
3 Trs^rv:::-  m moment  of  Inert in 
s Maas 
▼ Velocity 

N Spin  of  the  shell  about  its  longitudinal  axis 


b.  Geometric  Magnitudes  and  vectors: 

1,  j F-id  k - unit  vectors  pointing  along  -she  directions 
of  the  x,  y and  z axes  respectively. 

A the  unit  vector  having  the  direction  of  the  axis  of 
the  shell.  If  the  direction  cosines  of  A are  jf, 
a and  n,  A - jfi  + mj  + nk. 

5 the  angle  of  yaw  - the  angle  between  the  trajectory 
and  the  longitudinal  axis  of  the  shell. 

• 

v Ax  A t the  vector  Aguiar  velocity  of  the  axis  or 
the  shell. 


TO* 

4d 


=*2 


u /. 

X/* 


1/t 


c.  .Aerodynamic  forces:  Units 

. o 

V Drag  - the  force  acting  cm  the  projectile  along  the  mjf/t** 
trajectory  opposite  to  its  direction  of  motion. 

A A 

M Overturning  moment-  assumed  to  act  in  ths  plane  of  mjf/t 

the  yaw  and  assumed  to  be  proportional  to  sin  &. 

Thus,  M may  be  written  as  M ■ p sin  B,  where  p 
is  celled  the  moment  co-factor. 

O 

L Cross  wind  force  - a force  acting  perpendicularly  mjf/t 

to  the  trajectory,  in  the  plane  of  the  year: 
assumed  to  be  proportional  to  sin  &• 

A A 

Ifw  Tawing  moment  due  to  yawing  - when  a shell  yaws,  mjf  /t 

there  is  a torque  exerted  on  it,  the  axis  of 
which  coincides  with  the  axis  of  the  yawing  motion 
and  exerts  a moment  which  opposes  tha  yawing  motion? 

This  ma&cjiit  is  called  the  yawing  moment  due  to  yawing 
and  is  represented  by  Hw. 


rig 


Units 


K riaguue  force  - c orre spor ding  to  the  cross  vied, 
force,  when  a shall  spins,  there  la  a Magnus 
force  acting  on  it,  proportional  to  sin  6,  "but 
acting  perpendicularly  to  the  plana  of  the  yaw. 

J Magnus  moment  - the  moment  of  the  Magnus  force* 


•ARP  Rolling  moment  - the  moment  that  opposes  the 


XU- 

MA  VUW  PUIQX.A  ' 


f The  rolling  moment  damping  factor 

* AN 

* “ AN 


k The  cross  wind  force  damping  factor 


nrv  sin  o 

h The  yawing  moment  due  to  yawing  dangling  factor 
h « - H/B 

A The  Magnus  fores  damping  factor 


mv  R sin  6 


7 The  Magnus  moment  damping  factor 


AH  sin  5 


9 . 

These  damping  factors  all  hare  1/t  for  unit. 


^ . 


(*«<X  W3C 


Chapter  I 

Motion  with  Simplified  Force  System. 

To  simplify  the  problem  and  to  bring  out  the  essential  points  of  the 
theory,  we  shall  consider  first  a shall  subject  to  a force  system  con- 
sisting only  of  the  overturning  moment,  M,  the  drag,  D,  and  a spin  rertroy- 
Ing  couple  or  rolling  moment*  We  take 

tf  _ -J-  • 

r«  — p waU  C 

where  ^ ie  the  moment  co-factor  and  8 Is  the  angle  of  yaw.  We  assume  the 
magnitude  of  the  rolling  moment  to  be  AH  p where  A is  the  axial  moment  of 
inertia  of  the  shell  end  N the  spin-  The  axis  of  the  torque  AN  if  Is  that 
of  the  shell. 

Gravity,  the  cross  wind  force,  yawing  moment  due  to  yawing,  etc.,  are 

Iwa  vaso  ^IdOCUW  UUUbl<CU« 

Since  there  are  no  forces  transverse  to  the  trajectory,  the  trajectory 
is  linear.  Therefore  a coordinate  system  with  origin  at  the  center  of 
gravity  of  the  shell,  with  its  x axis  pointing  along  the  trajectory,  with 
its  y axis  In  a vertical  plane,  and  the  z euls  pointing  to  the  right,  while 
not  Galilean  with  respect  to  translation,  is  Galilean  eo  far  as  angular 
motions  are  concerned.  See  Figure  1. 


y 


The  unit  vectors  having  the  directions  of  the  x,  y,  z axes  are  designated 
as  i,  j,  k,  respectively,  (These  vectors  as  well  as  all  other  vectors 
ere  written  with  heavy  letters.)  we  define  the  relevant  quantities. 

B B = moment  of  inertia  about  a transverse  axis  through  the  center  of 

gravity. 

A / A - unit  vector  in  the  direction  of  the  axis  of  the  shell. 

w w - the  resultant  vector  angular  velocity  of  the  longitudinal  axis 

about  a transverse  axis  through  the  center  of  gravity  (also 
called  tha  cross  spin). 

tt  H * total  vector  angular  momentum  of  shell  with  respect  to  0. 


By  reason  of  ths  symmetry  of  revolution  vhich  the  shell  is  assumed  to 
possess,  the  total  angular  racaantura  is  the  sum  of  tvo  components: 

(a)  ths  component  about  the  axis,  vhich  has  the  direction  A and 
is  represented  "by  the  vector  AKA.,  and 

(b)  the  component  transverse  to  the  axis,  vhich  is  represented 
by  the  vector  Bw. 

M. 

x bi"»  vector  cross  spin  v is  a quad  in  magnitude  to  A and  is  perpendicular  to 
A and  A . He  thus  have 

V > A x A .. 


The  total  angular  momentum, 

E - AHA  + B [~A 

Thr  vector  equation  of  motion  is  that  the  rate  of  change  of  angular 
G momentum,  H,  Is  equal  to  the  Impressed  vector  couple,  designated  by  Q. 

If  ve  take  N as  variable  ve  fixd  that 


[A,  A] 


r ■ .. 

+ B |^A  X A 


(1) 


H - AKA  + AHA  + B 

• • • 

The  A x A term  in  2e.r0  since  A and  A are  parellel.  Thus 

E = AKA  + AHA  + B I A X A I w G 


]• 


• •"1 
A X A |- 
- J 


On  the  present  assumption,  the  vector  couple  G consists  of  tvo  terms,  the 
couple  M end  rolling  assent  ANf"*  . Let  us  see  how  we  nay  express  them 
as  vectors. 

M is  equal  in  magnitude  to  4 sin  & and  the  axis  of  the  couple  is 
obviously  perpendicular  to  the  trajectory  and  hence  to  the  unit  vector  i, 
and  to  the  axis  of  the  shell. 

r~  -i 

Consider  the  vector,  [i  z Aj  i this  vector  is,  by  definition  of  the 
vector  product,  perpendicular  to  1 end  A and  has  the  magnitude  sin  5. 

Hence,  the  moment  M is  expressed  aa  the  vector,  4 jjL  x a]  • The  rolling 

moment,  AE  T , has  by  definition  A as  its  axis  and  the  moment  is  assumed 
to  be  negative.  Bence  ve  represent  it  vectorially  as  -AN  p A. 

Replacing  ® in  the  preceding  equation  by  itt  tvo  ccaqponents,  uTi  x A I 
and  -AH  p A,  the  equation  becomes  *—  -* 

(la)  H = AHA  + AHA  + B | A x A~|  - \x  fi  x A~1  -ANf  A ; 

1 — — ^ L-  -J 

this  Is  the  vector  equation  of  motion. 

Differentiation  with  respect  to  the  time,  t,  is  indicated  by  the  superposed  dot. 


2 2 


**In  ballistic  notation,  ^K^pd  v 


It  should  he  noted  that  since  A is  a unit  vector,  A 1b  necessarily  per- 
pendicular to  A,  while  [a  x X]  and  jji  x aH  are  perpendicular  to  A by  the 
definition  of  a vector  product.  Hence,  if  we  take  tha  scalar  product  of  A 
into  equation  (l),  all  the  term  vanish  except  AN (A.  A)  and  -ANT  From 

this 

AN  - -AN P 

and 

ANA  - -AN  P A. 


By  virtue  of  this,  it  is  obvious  that  (la)  may  be  written 


(1.01) 


ANA  + B 


J^A  x A J « p jjL  x Aj 


If  the  direction  cosines  of  A with  respect  to  the  x.  y,  z axes  are  f,  m,  and  n, 

A ■ jfi  + mj  + nk 

A = jfi  + mj  + nk 
»• 

A ■ fi  + mj  + nk 

and 

A x A J » (jfi  + mj  + nk)  x ( jfi  + mj  + nk)  . 

Upon  performing  the  vector  multiplication,  remembering  that 
i x 1 ■ 0,  i x ,1  - k,  J x 1 > -k,  etc.,  ws  have 
A x A * (mn  - n£)  i + (nf  - jfn)  J + (jfm  - ajf  )k  . 

For  i x A,  we  have 

i x (jfi  + mjJ  + nk)  ■ - nj  ink  . 

EqUouIou  (1.0%)  written  out  in  full  becomes 


(21.01X) 


AN(fl  + mj  + nk)  + B j (inn  - nm)i  + (njf  - jfn)J  + 

+ (fm  - mjf)k  j - p j~-  nj  + mk  J . 

In  this  theory,  the  assumption  is  made  that  the  yaw  is  small.  If  the 
yaw  is  small,  then  the  i component,  the  magnitude  of  which  is  cos  d,  is 
slvsys  taken  as  unity,  and  therefore  the  i component  is  not  involved  la 
this  theory. 

Since  the  i,  J and  k components  are  independent,  there  are  three 
equations,  one  for  each  component.  We  take  the  4 end  k components : 

AHa  + B (nj  - jfn)  - - pn, 

• • 

AXn  4-  B (jfm  - mjf)  * + pm  . 


(2) 

43) 


(j  component ) 
(k  component) 


8 **  • 2 •* 

f is  equal  in  magnitude  to  coa  6 « 1 - and  ( m - 6 -68# 

•* 

Nov,  m and  n are  of  the  order  of  6.  Bence,  in  neglecting  such  terms  as  n /, 
V*;  are  neglecting  a term  of  the  order  of  (66^  - 6~  8)  in  comparison  vith  6 
and  in  taking  i - 1,  ve  are  neglecting  a term  of  the  order  of  52  in  cocrpar- 
iaon  with  unity.  • We  nov  t&ke  ( - 1 and  neglect,  the  terms  in  jf  . 

Equations  (2)  and  (3)  then  become 
(2.1)  AHm  - Bn  « - 

(3«l)  ANn  + B»  ■ * jin  * 

Multij  ‘ y (2.1)  by  i ■ -1  and  then  subtract  from  (3*1)  vith  the 

result, 

AN(-im  + n)  + B(a  + in)  » ji  (k  ♦ in). 

Since  -(i)2  - -l(-l)  - +1,  the  first  term  may  be  written  -lAH(m  + in) 
and  the  equation  becomes 

(4)  -iAlf(m  + in)  + B(S  + in)  - u (a  + in). 

i)  We  nov  substitute  tj  for  (m  + in)  and  call  r,  the  complex  yaw. 

The  complex  yaw  tj  Is  represented  on  a plane  perpendicular  to  the 
trajectory  as  indicated  in  Figure  2. 


GxiS  of 

Figure  2 


The  axis  of  reals  points  upward  in  a vertical  plane  along  the  y axis 
and  the  axis  of  imaginarles  points  to  the  right  along  the  z axis. 

With  tj  substituted  for  m + in,  (4)  becomes,  upon  division  by  B 

(4.1)  tj*  - i Q tj  - | ij  - 0, 

Q if  Q is  used  to  represent  AN,/!}. 


4 


To  solve  this  equation,  it  is  convenient  first  to  eliminate  the  f| 
tern.  Let  us  write  it  for  future  convenience  in  the  form 

(4.2)  rf  - 3LEt|  - Ftj  « 0 . 

Ve  new  make  the  substitution 


n - y exp 


si  ( Edt 
o 


Upon  differentiating  once,  we 

■*=»  - — j 

have 

f 

q » y exp 

t 

% i ( Edt 

d 

y dt  ' 

[ 

|e*3> 

*»  ^ 

j Edt 

D 

o 

— — 

> 

y e*P 


1 • / 
2 1 


1 


\ Lit:  | + y exp 

1 J 


i J .1 

f i ) Edt  (f  i E) 

o 


am 


1 U 

1 i f Edt 

2 i 


1 


y + | i E y 


Hie  second  differentiation  given: 

r 

..  a j 

* - at  ) exp 


+ 

1 1$  Edt 

(y  + | iEy)  + exp 

| i ; Edt 

o 

1 

>T 

o 

exp 


exp 


i 

I1 


Edt 


_ i 

i if  Edt 


| i E{y  | iEy)  + ? + ^ + . 

y +.  (§  IE  + | iE)jr  + (-  J E2  + 4^)yl  . 


Unon  substituting  for  q,  q>  and  V *n  (^.2)  it  becomes 


exp 


1 t 
iij  Edt 


- exp 

|i(Edt 

« 

o 

1 ll.+  ^ + ^ + ^)J0  “ + I iESy) 

/E^  _ 1 | _ 


- Fvl 

J 


and  the  equation  for  y is 

(fc*3)  y + (tj*  - p + \ i E)y  - o . 

5 


* 


rwuho*. 


f 


Upon  substituting  for  js  and  F their  values  as  obtained  frpm  (4.1) 
and  for  the  present  taking  ft  as  constant,  that  is,  neglecting  E in  (4.5), 
we  have 

(4.4)  y + <tj-  - |)y  - 0 . 

This  may  be  written 


if 


aS2 

4Bp 


fl  - Jin 

• 

2 

. 1 

L 

%r  - * 

'S' 

A " a 

u 


f'b  c\ 


designating  ijl  - l/s  by  a,  we  finally  get 

••  /H  \P 
y + VjjQ/  y « u . 


For  convenience,  following  Kelley  and  McShane  , we  rewrite  (4.5)  as 

(4.6)  y - q2y  « y + (iq)“y  = <f* 

and  make  another  substitution 


From  this 


r * 


L22ZJL 

dt 


V 

*■  . 
y 


-•>  7 


- 


ry  + 


2 

r y 


and  (4--. 6)  becomes 


. 2 2 
ry  + r y - q y 


0,  from  which 


(4.7) 


If  f is  small  in  comparison  with  r,  we  should  have  two  solutions 
for  r in  this  equation,  one  nearly  equal  to  +q,  and  the  other  nearly 
equal  to  -q.  Hence  let 


r = + 5 t c • 


By  substituting  in  (4.7),  it  appears  that 


2 


+ q + C + q‘-  + 2qe  + €**  - q ® 0, 
from  which,  if  v*  neglect  k and  c , 


. A 

2q  * 


d r e + q - *;5£ 


^Kelley-hcShane  “ ERL  Report  446,  1944T 


**Kelley-McShsne  refer  to j H.  Jeffreys,  Proc.  loadon  Math.  Soc.  (2),  vol.  23 
(1923)  p.428,  aud  G.  Wentzel,  Zeltscbrift  fur  Physik.  vol. 38  (1926)  p.  518. 
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From  the  definition  of  r, 

rdt  * d( log  y) 


f rdt  * log  y - log  y , 

6 

t 

From  this,  log  y ■ log  y^  + ( rdt 

6 


y - yQ  « 


-ye 


£ (±q  - ^q)dt 


•imus,  the  two  solutions  for  +q  end  -q,  respectively,  ere 

t . 

y - yx  exp  ) (q  - -^)dt 


y « y2  exp  £ (-q  - ^)di  • Thus,  y^  and  yg  are  constants. 


From  (4.5)  cuid  (4.6)  we  have 

2 ,iaj,a 
5 = (-y) 

Ifle 

+q  - + -g- 


ifle 

-q  « - -5-  • 


Furthermore 


2q  " 2±a  a 

Since  A and  B are  constant, 


± (2.  . 

2 'Q  + o 


e • 

a n 


^ which  we  designate  by  - ( , 

where  jf  is  called  the  rolling  moment  damping  factor*  We  have  then. 


% - 


substituting  for  q and  q/2q, 


for  the  two  solutions 


y = yx 


y j'o  exP 


5 f (iflo  + / 


• !>  *] 

|§  j(-  mo  + ( - 2)  atj  . 


(5) 


Recalling  that  r = y exp  | % i \ ftdt  | and  that  the  general  solution 
is  the  sum  of  the  two  solutions,  we  have 

P1  P2 
Tl  » e + ^ e 


where 


and 


si  ■ H f 


ifl(l  + o)  + jf  - 


dt 


p2  = 1 J jiot1  - O)  ♦ if  - 1 1 at  . 


In  (5) , K,  and  are  arbitrary  complex  constants  determined  by  the 
initial  conditions. 


Interpretation  of  Solution 


For  those  who  have  not  had  experience  in  the  field  of  complex  quantities, 
an  interpretation  of  the  solution  (5)  is  offered. 

Consider  the  series  expansion  of  ey. 


y , y 

e » 1 + y + 7jj- 


+ *? 


+ etc. 


Upon  substituting  lx.  for  y in  this,  we  have 

k 


ix 


1 A i V 


X2  ix5  X _ _ , 

x = 2f  " ~$T  + WT  ^ w “••• 


ix^  x^ 


2 ^ 6 

27  + 77  - 57  + 


+ 3 


x5^*5, 
x * jT"  5T  + *** 


ix 

The  first  term  on  the  expansion  of  e is  the  series  for  cos  x while  the 

8 


expression  in  the  second  bracket,  is  the  series  for  sin  x.  Hence  vs  have 

ix 

e «*  cos  x + i sin  x • 

For  the  case  when  p Jj.  x A~j is  the  only  component  of  the  force  system 

and  & m ( m o f the  solution,  as  ve  have  seen,  is  of  the  form 

|iO(l  + o)t  |ia(i  - a)t 

tj  =>  K^e  -t  Xg  e • 

Lee  us  write  this  as 


n - 


and  let  us  also  assume  for  simplicity  that  and  are  both  real. 


Since  a is  the  real  part  of  r;  and  n the  imaginary  part  of  i),  then  in 

■j  y 

view  of  the  interpretation  of  e , (6)  may  be  written  a a? 
f»°h  cos  ®jt  + Kg  cos  o^t 


•^n  = Kj^  sin  ®jt  + Kg  sin  «gt 


consider  an  epicycllc  motion  of  the  following  character.  Let  there  be 
a circle  of  radius  It  as  shown  in  Figure  J »n«*  let  a point  P move  cm  it  with 
a constant  angular  velocity  a^.  With  P as  a moving  center,  let  there  be 
another  circle  of  radius  JC,  and  a point,  Q,  moving  on  it  with  a constant 

jO  - o 

angular  velocity  axj . v ~ u>,t 


/ 


'»  2 


Figure  3 


If  fl  in  o is  the  axis  of  m's  and  9 » jc/2  is  the  axis  of  n's,  we  have 
for  m and  n at  any  time  t, 

3 *a  CCS  OC^t  + Xg  COS  Oigt 

n ■ sin  c&jt  + Kg  sic  a^t  . 

But  this  is  the  same  as  given  hy  (6*1)  and  the  result  establishes  the 
cpicyclic  Cwu;<»:l«r  ui  tue  motion  given  in  \t>)* 

In  cases  of  practical  interest,  neither  a nor  f is  zero?  Let  u® 
rewrite  P,  and  P0  of  (5)  as 

1 c. 

P1  - i | J 0 (1  + a)dt  + | ( (f  - a/o)  dt 

P2  “ 1 I j “ (i " °)dt  + i ( ” ®/<0  dt  • 

The  first  parts  of  P^  and  Pg  are  imaginary  and  the  second  parts  are 

real.  In  view  of  the  fact  that  o is  no  longer  taken  as  a constant,  the 

angular  velocities  and  au  are  no  longer  constant.  The  real  second 
parts  indicate  that  the  amplitudes  are  no  longer  either. 


So,  in  place  of  (6.1),  we  write 


c rr  r — t 

^ K^cos  jjj  ^n(l+o)dtJ  + KgCos 

“ r ~n 

I j 0(l-o)dt 

> 

J” 

r.  r -J 

7J  1 (jf-&/o)dt 

/it,  sin 

”fn(l+a)dt  | r K_sin 

L?*  J * 

1 J a(i-e)dt 

}- 

jj  J 
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Chapter  II 


Equations  of  Motion  with  a Spinning  Coordinate  System* 

The  Spins  Caused  by  the  Cross  Wind  Force,  the  Magnus  Force  and  Gravity* 

The  Spinning  Coordinate  System 

In  the  preceding,  ve  have  followed  the  methods  of  Fowler,  Gallop, 

Lock,  and.  Richmond  ia  deriving  equations  (4.1)  and  (4*5)  and  the  methods 

#.i  i j v—  ir»n  ...  A *Ua  aa1..4.4  a..  t«  j A.  J 

A avwvw.  wj  ■ - *-  • * • j * mm » »ww« » '«/  *» v* ««* *-*Q  »/  wamvawm * xm  »4> 

the  motion  In  the  presence  of  a cross  wind  force,  and  other  forces  trans- 
verse to  the  trajectory,  we  depart  from  Fowler,  Gallop,  Lock  and  Richmond, 
and  adopt  a coordinate  system  which  is  attached  to  the  actual  trajectory* 

Art  In  Chapter  I,  the  unit  vector,  i,  points  along  the  trajectory,  the  unit 
vector,  J,  is  initially  in  a vertical  plane  through  the  trajectory,  and 
the  unit  vector,  k,  Is  initially  horizontal.* 

Aa  a result  of  the  cross  wind  force  and  other  forces  having  components 
transverse  to  the  trajectory,  our  coordinate  system  will  have  a component 
of  spin,  to,  about  some  axis  transverse  to  the  trajectory  hut  no  spin  about 
the  trajectory.  By  virtue  of  the  spin,  the  coordinate  system  will  not  be 
a Galilean  one.  kTe  proceed  to  derive  first  ra  and  then  uhe  equations  of 
motion  in  this  spinning  non-Galilean  coordinate  system* 

Consider  a body  of  mass,  m,  moving  with  velocity,  v,  in  the  direction 
indicated  by  the  unit  vector,  T»  Its  vector  mooentum  is  evidently  mvr. 

The  rate  of- change  of  momentum  or  the  vector  force  is  made  up  of  two 
components,  ran  and  snrir.  The  first  component,  parallel  to  v,  is  the  sum 
of  the  drag  and  the  tnmgential  component  of  the  force  of  gravity,  while 
the  second  component  ie  a force  perpendicular  to  t*  Inis  1b  apparent 
since  r Is  a unit  vector  and  the  only  possible  change  in  t is  a change 
In  direction.  See  Figure  4. 


Figure  4 

If  we  designate  the  vector  force  perpendicular  to  r by  F , we  have 

T ^ 

T > — — - • 

_________________________  nxv 

* 

The  deviation  of  k from  the  horizontal  direction  and  J from  the  vertical 
plane  will  be  discus  sod*. in]  .CJiapber  % It  will  be  shown  that  theee  deviations 
are  small. 


is  equsl  in  magnitude  to 


Thus,  the  angular  velocity  of  the  trajectory,  as, 


F 

a 

mv 


As  a vector,  it. must  be  perpendicular  to  T^and  also  r.  Hence,  the 

spin  of  our  coordinate  system,  (d**txt«*t2c--^. 

rav 

TP. 

**  D p 

It  is  apparent  that  — — may  be  replaced  by  — . where  F is  the 

mv  mv  7 

resultant  force  on  the  shell,  since  the  tangential  cnrrmnnent  of  F ashes 
no  contribution  to  the  vector  product. 

In  the  future,  vs  shall  replace  r by  1 since  by  definition,  1 is  the 
unit  vector  pointing  along  the  trajectory.  So  that 


CD  - i X $. 


i X 


mv 


Fro®  page  2,  we  have  as  the  equation  of  motion 


«•  _*  I ..i 

Am  + Aim  + a (_A  x A_|*  0 • 

The  superscript  . indicates  time  derivatives  in  a Galilean  coordinate 
system.  If  superscript  primes  indicate  time  derivatives  in  our  chosen 

coordinate  system,  the  Job  is  to  express  N,  A and  A in  terms  of  N*  . 

* • a f a ... 

IX  y A eUlU  U)» 

As  the  first  step,  we  state  the  following  theorem: 


Figure  5 

In  this,  rj  is  the  rate  of  change  of  the  vector,  r,  as  measured  in 

space  1 (non-Galilean)  with  axes  x1,  y^,  Suppose  this  space  is 

spinning  with  respect  to  Galilean  space,  2,  with  an  angular  velocity,  cd, 

aa  indicated.  It  is  apparent  that  the  velocity  in  space  2 of  the  point 
P at  the  end  o?  r caused  by  the  spin  o>,  will  be  a x r.  It  may  be 


is  to  be 


rigorously  shewn  that  tbs  rate  of  change  of  r in  space  2,  , 

obtained  by  adding  this  term  to  . Hence,  sinee  space  2 is  Galilean, 


rg  ■■  r£  + o)  x . 


jjYom  this  theorem,  it  follows  directly  that  A ■ A'  + © x A and  by  a 
second  application  that 

A"  + a x A*  + os1  xA  + fflXA'  + as  x jjn  x A~j  » 


‘A* 


■ A"  + 2 © x A'  + ©*  xA  + «x[»xa] 


Although  © has  no  component  about  the  trajectory,  it  does  have  a 
component  os  sin  6 about  the  axis  of  the  shell,  which  makes  an  angle  b 
*rith  the  tangent  to  the  trajectory*  We  shall  neglect  © sin  6 in  compar- 
ison iflth  N.  -~S<5  wS  take  "the  spIA  N unchanged  in  our  coordinate  system. 

Upon  making  the  substitutions  indicated  above  in  (l),  vs  obtain 
AN 'A  + ANA'  + AN©  x A 


Co\ 
f~  J 


+ B 


|A  x A"  + Ja  x [«'  x A]J  + 2 j~A  x x A'^Jj 

«tQ  j * ° 


+ [a  X |~C5  X jjB  X A 

To  simplify  this  result,  we  make  use  of  a theorem  of  vector  analysis, 

jTx  [1  x Cj  j = B (A«C)  - C (AtB)  . 

For  the  two  triple  and  one  quadruple  products  indicated  by  the 
superscripts  I,  II  and  III,  respectively,  we  obtain 

(I)  !"a  x Jos'  x a”!  * ©'(Aia)  - A(a?«d')  « »'  - (A^ffl')A 

U.  -4 

(II)  2pA  x F©  x A'!]  « 2tt>  (a*A')  - 2A'  (A*®) 

L ^ ^4 J 

(III)  [a  x[  © X o* * aj11  - uhA'  | © X *D-|wi  (a-®) 


In  treating  the  term  ( A*®'  )A  of  I,  we  consider  two  cases.  The  first 
case  is  the  one  where  © arises  frc/m  aerodynamic  forces  and  is  therefore 
proportional  to  $.  Since  in  the  final  A,  we  consider  only  the  J and  k 
components  which  are  themselves  proportional  to  5,  It  is  apparent  that, 
^The  use  of  the  superscript  dot  ie  appropriate^  since  space  2 is  Galilean. 
- "The  name  Coriolis  is  associated  with  a similar  transformation*  See 
"Classical  Mechanics",  pp.  136-137>  “by  Goldstein. 

1? 


(1*03) 


for  this  case,  (a*®5  )a  is  at  least  of  order  5&J  and  may  be  neglected  in 
this  first  order  theory * In  the  second  case,  © arises  from  gravity  and 
© and  ©3  are  always  perpendicular  to  the  particle  trajectory;  to  the* 
actual  trajectory,  they  will  be  so  nearly  perpendicular  that  the  product 
(A*©*)  will  be  of  the  order  8 and  (a*©' )a  of  the  order  of  8^.  This  also 
is  neglected. 

In  II,  since  A is  a unit  vector.  A'  is  necessarily  perpendicular  to 
A*  Bence,  the  scalar  prodxict  (A*A')  = 0.  The  spin  a is  perpendicular 
to  i,  hence,  the  product  (A*©)  is  at  least  of  the  order  8 and  A'  is  of 
the  ordftr  8’.  Avrirpsainn  -Pa i / a »/rA  ts  at  least  of  the  order  65*  snd 

is  thus  omitted  in  this  small  yaw  theory. 

, ,,  , . . - r-  -I  . r~  *1 

in  in,  one  scalar  produce  ^a*  |©  x a l ; * u,  since  |®  x A I is  per- 
pendicular to  A.  *“  “* 

In  view  of  these  considerations,  (1.02)  may  be  written 


If  the  spin  © is  caused  by  aerodynamic  forces,  © will  be  proportional  to  8 

E~1  . 2 

t x A j (A*©0  will  be  at  least  of  the  order  8 and  thus  may 

oe  uegi.ee  bea. 


We  proceed  to  compute  the  spins  , ©.. , 

R ' 

wind  force,  the  Magnus  force  and  gravity, 

the  resultant  aril  n m m m r ® + S . 

i\  a g 


cx  , and  © caused  by  the  cross 
a-  S 

respectively,  thus  obtaining 


The  Cross  Wind  Force 


The  cross  wind  force,  L , is  defined  to  be  a force  which  acts  per- 
pendicularly to  the  trajectory  in  the  plane  of  the  yaw.  For  this  treat- 
ment, ills  further  assumed  to  be  proportional  to  sin  5.  So  we  take 

L « A sin  8, 

and  call  A the  cross  wind  force  co-factor. 


Consider  now  the  vector  difference 

(A  - cos  8 i);  (see  Figure  6). 


*In  ballistics,  we  write 

2 2 

L p p d u sin  8. 


(See  Hayes1 
Ik 


Elements  of  Ordnance^  p.  Ul2) 


It  is  obvious  that  (A  - cos  5 l)  has  the  magnitude  sin  b and  that  it  is 
perpendicular  to  the  trajectory,  the  direction  o.f  which  is  i.  Accordingly, 
a correct  vector  representation  of  the  cross  wind  force  is 


>.  (A  - cos  5*1)  . 


Frciu  (?)  it  appears  that  the  spin  © . of  the  coordinate  system, 
caused  by  the  cross  wind  force,  is 


4 „ * (A  - cos  5 1)  J * 

a ww 

mv 


j\_ 

mv 


L‘xAJ 


mv 

T 


mv  ^ 
sin  5 QH 

■■■  *M| 


* a a • o jji  a i | * = 


.“5  l 

LJ? 


Upon  replacing  — or 5 — t by  k,  we  have  for  the  spin  © caused 

* 0 mv  mv  ein  % ’ x 

by  the  cross  wind  force, 


©. 


Magnus  Force 


Hie  Magnus  force  is  a force  which  arises  from  the  interaction  of  the 
boundary  layer  of  a spinning  shell  and  the  wind  stream  when  the  shell  has 
an  angle  of  yaw.  Consider  a tennis  ball  or  a baseball  under  the  conditions 
shown  in  Figure  7*  As  a result  of  the  interaction  between  the  wind  stream 
and  the  boundary  layer,  the  velocity  at  the  top  surface  of  the  bail  will  be 


•wind  velocity 


less  than  the  velocity  on  the  bottom  surface  of  the  bell.  Associated  with 
this  will  be  a higher  pressure  over  the  top  thus  there  is  on  the  bottom 
producing  a force  ifaich  accelerates  the  ball  in  a downward  direction.  This 
is  what  happens  when  a tennis  ball  is  given  a top  spin. 

Consider  a shell  spinning  as  shown  in  Figure  8.  As  a result  of  the 


Spit?  of  9h*l/ 


Vu 


SC 


’ tfofm ’£  fe^ce 


Figure  8 
15 


pressure  distribution  mentioned,  there  will  be  a Magnus  force,  K,  acting 
perpendicular  to  the  plane  of  the  yaw,  as  indicated,  which  should  be 
approximately  proportional  to  the  spin,  K,  the  velocity  of  the  shell,  v, 
and  to  the  sine  of  the  angle  of  yaw,  5.  If  the  factor  of  proportionality 
is  represented  by  fm,  we  may  write 

K K*  = fmvN  sin  6 (in  magnitude) 

and  by 

I — — l 

ftavN  I A x i | as  a vector. 

X sfe  shall,  represent  fN  in  the  above  by  X and  call  X** ***  the  Magnus  force 

<i«mping  factor.  The  contribution  to  i caused  by  K i«  therefore 

i = Ia xj)  , w i~~i  . 

mv  ti J 


The  spin  ea^  caused  by  the  Magnus  force  is 

j 1 a 1 i — a j i x |^A  x i | | . 


Gravity 


g The  force  of  gravity  is  numerically  equal  to  mg  where  g is  the 

acceleration  caused  by  gravity.  This  force  has  a component  perpendicular 
0 to  the  trajectory  equal  in  magnitude  to  mg  cos  0,  where  0 is  the  inclina- 
tion of  the  trajectory  to  the  horizontal.  Therefore,  in  our  coordinate 
system,  the  component  perpendicular  to  the  trajectory  will  be  represented 
by 


- mg  cos 


It  follows  from  (7)  that  the  spin  u>_, 

o 


i x 


(-g  cos  e j) 
v 


6 1 


since  0‘ 


- mg  cos  0 
v 


caused  by  gravity  , will  be 

k 


(8) 


The  Resultant  Spin  of  the  Coordinate  System 

The  spins  caused  by  the  cross  wind  force,  the  Magnus  force  and  gravity 

uac  respectively , 

f °>k  - * L1  * a3  * 


< 


"X 


a> 


« X i x 


0’  k 


If11]  - 


i mV ° 2 

*In  oallis*bic  notation,  w*  write,  Magnus  force»K«K^pd  sin  6(in  magnitude)* 


* y!ne  damping  factor  X is  to  be  distinguished  from  the  crosswind  force  co- 

factor X* 

***In  view  of  the  small  rotation  of  cur  system  about  the  trajectory,  the 
component  of  the  gravity  force  will  not  be  exactly  represented  by  - mg  cos  0J 
There  will,  in  general,  be  a sms  13.  k component. 
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Chapter  III 


The  Yawing  Moment  Due  to  Yawing. 

The  Magnus  Moment. 

The  Equation  for  the  Complex  Yaw. 


We  have  so  far  considered  the  following  elements  of  the  force  system: 

Drag 

Overturning  Couple 

AvI  a l ^ 

Cross  Wind  Force 
Magnus  Force 
Gravity 

The  last  three  forces  contribute  to  the  spin  of  the  trajectory  and, 
hence,  to  the  spin  of  our  coordinate  system.  In  this  chapter,  we  consider 
the  yavlng  moment  due  to  yawing  and  the  Magnus  moment,  and  derive  the 
equation  f >r  the  complex  yaw. 

Yawing  moment  due  to  yawing 

# 

H The  yawing  moment  due  to  yawing,  designated  by  Hw  , is  the  moment 

which  opposes  the  angular  velocity  of  the  axis  of  the  shell. 

h We  define  a yawing  moment  due  to  yawing  damping  factor,  h,  by  writing 

ftr  » - h B v • 

As  we  have  seen  in  Chapter  I (p'\gr?  2),  the  angular  velocity  of  the  axis  of 

the  shell  w may  be  expressed  as  (a  x A).  Thus,  if  we  are  to  take  account 
of  this  moment,  we  include  on  the  right  hand  of  equation  (1.03)  on  page  14, 

- h B (A  x A) 

in  addition  to  the  overturning  moment,  p |l  x A I , and  the  axial  couple, 

-AN  PA.  U J 

• r 

As  has  boon  mentioned,  A is  the  angular  velocity  referred  to  Galilean 
axes,  where as,  in  this  problem,  we  refer  the  motion  to  axes  turning  with 
the  trajectory.  We  thus  have  to  express  A in  terms  of  A*  and  thv  reoultsct 
angular  velocity  of  the  trajectory,  es  As  we  have  seen  in 
equal  to  A*  + <s  x A,  and  thus 

A x A ■ A x A'  + A x [_•  x A^  « 

As  we  have  also  seen  in  Chapter  II,  the  triple  product  A x 
expressed  as  to  - (o  s A)a«-  Hence 

-hB  [a  x Al-  -hB^Ci  x aL  + co  ~ (m  • A)A 
*In  ballistic  notation,  Hw  Is  expressed  by 

Kg  p J*  vw 

in  which  K^  is  the  yawing  moment  due  to  yawing  coefficient 

vector  angular  velocity  of  the  axis  of  the  shell. 
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Chapter  II,  A is 


, w x A ; may  be 

4=  =4 


| A 

L 

} 


and  w is  the 


The  term  (®  * A)A  will  "be  omitted  since  it  is  of  the  order  5 in  the 
4 and  k components.  Thus  we  have,  including  the  moments  sc  far  considered, 
for  the  right  hand  side  of  (1.03), 


|7  x.  aJ  - AHpA  - hB  ^|^A  x A'J  + <o  J 


On  page  15,  it  appeared  that  the  vector  Magnus  force  is  proportional 
to  A x i.  The  moment  of  this  force  Is  perpendicular  to  the  force  am  to 
the  axis  of  the  shell  and  is  the re foie  properly  £iven  as  proportional  to 

] — — i 

A x j A x i j . We  assume  a proportionality  factor  of  Aliy,  where  / 1b  the 
Magnus  moment  damping  factor,  and  write  the  Magnus  moment, 


J*  « AH7  |X  x [a  x i^j~ | . 


This  term  is  to  be  added  to  the  right  hand  of  equation  (1.03)  in  addition 

to  the  moments  mentioned  above. 

While  it  is  commonly  assumed  in  the  small,  yaw  theory  that  the  Magnus 
force  and  moment  are  strictly  proportional  to  the  spin,  there  is  little 
experimental  eviclence  to  back  up  this  assui^ption.  What  evidence  there  is 
indicates  that  the  force  and  moment  are  not  exactly  proportional  to  the  spin. 
On  the  other  hand,  the  evidence  indicates  that  the  Magnus  force  and  moment 
are  proportional  to  8 for  small  yaws. 

Including  all  the  moments  so  far  considered,  the  right  hand  side  of 
(1.03)  is 

- D-0  - AN  Pa  - hB  ^ |^A  x A^J  + &J  + ANy  J^A  x |~A  x . 

If  the  term  - AN  pA  is  cancelled  against  AN 'A,  (l,03),  see  page  14, 
becomes  v 

f—~  — t _ f h—  — ^ r~  — t / 

(1,04)  AN\A'  + j CD  x A j)  + j A x A”j*  - (A  * ®)  [_®  x Ajy  = 

- p Jjl  x A^j-  hB  ^ A x A' J + ®^>  + AN 7 jji  x jji  x J 
with  the  spin,  as  given  in  equation  (7)  of  Chapter  II  by 


£jl  x A "|+  A ix|^.x  + 0 * k 


We  are  to  expand  (l*04)  in  terms  of  the  4 so*  k components,  the  compon- 
ents of  the  complex  yaw,  and  omit  the  i components , if  any,  since  in  this 
small  yaw  theory,  jf  - 1. 

*1”  ballistic  notation,  we  write  J = Kjp  N v sin  8 (in  magnitude)* 
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Before  carrying  out  this  expansion,  it  is  convenient  to  expand  the 
vectors  and  the  vector  products  in  (l«04). 

For  m,  on  carrying  out  the  indicated  vector  multiplications,  ws  find 
with  . 

A - jf  1 + mj  + nk, 

55  i=  (=an  -t  Xs)J  + (an  + Xn  + Q*)k  « 

If  the  J and  k components  of  the  vectors  and  vector  products  In  (l*f)4) 
are  indicated  by  the  notation  ( ) , ».  we  have 


' 0 0 


,£’)»  v - m*  J + n •£. 

» j ~ 


I ! 


l 

X 


(#iA),  . ■ (am  + Xn  + 9 ' ).J  + (an  - >«)k  . 

J > E 

With  the  usual  approximations. 

(k  sr  . m -n»»  .1  + a"  k 

'J,  k ~ - 

(®*)j  ^ - (-an’  + Xjn')J  + (am1  + Xu'  + 3t5)k  . 

In  treating  the  term  (a  * A)  |jp  x a”]  we  need  to  Include  in  m>  only  the 

nentt  e’k  since  « will  be  proportional  to  B for  the  other  components  and  the 

2 

term  itself  will  be  of  the  order  6 for  these  components,,,  .. 


With  « - e% 

5 


f(*v  a>[va1}j,  k 


nd,cJ 


On  the  right  hand  side  of  (l.Oh); 

fa*  All 


(A  CA  1 Ojj,  * ' * 


mk  , 

+ rs^k  , 


nk  c 


/ i Ak\ 


Using  these  egressions  for  the  vectors  and  vector  products  in  (i.Os), 

AM 

dividing  both  rides  by  B,  expressing  the  quotient  -yr  by  fi,  and  collecting 
the  terms  in.  and  k,  respectively,  we  obtain 

O LL 

S(m*  +am  + X?:.  + 3’)  *>  n"  + (-an’  + Xm’)  - n3*  « - £ n + ha*  + ban  - hXott-Q7BU 

O 

f 

fl( h'+an  **  >*i)+  mM  * (am*  + Xu*  * 3*)  » » - hm’-bam  - h3\n  - h3’  + Oyn  v 


... 


We  isultiply  the  first  of  these  equations  by  -i  --  - "f--l  end  ade  to  the 
second  with  the  result: 

fi  |j-  im'  + n'  - k (iin  - n)  + X(-m  - in)  » 1 6 8 | 

•r  m,c  + in 90  r x(m5  + in0)  - *(im’  - ns)  + in0#2  + 0M  = 

■ ^ (m  + in)  - h(m’  + in’)-  liic(m  + in)  + h 7v( ±m  - n)  - h6-  - Ciy(im  - n)  . 
Replacing  (m  + in)  by  tj,  the  complex  yaw,  03  in  Chapter  I,  and  arranging 

olic:  t-rjinus , <rc  have 

TJrc  + »CT| 9 - i?M:  - iwT}S  - iilKf; 

-AXt]  - ^ q + hri*  + h*ci^  - ih>vT^  + iOyTj  ■ iQ0!  - 0”  - in 0’2  - hQ° 


4l-p6)  tj”  “ i(ft  + Ik  + ih  + X)tj 


-tF 


+ iftie  + flX  •>  ift7  - hK  + 


iO0°  - e"  - in  0°  - h0J  « 


ih>rji 


Chapter  IV 


The  Solutions  of  the  Equation  for  the  CgctdIsx  Yaw.  Stability. 

In  this  chapter,  ve  derive  the  solutions  of  the  equation  for  the  complex 
yaw  (l»06).  we  first  consider  a particular  solution  of  the  equation*  We 
insert  the  partioniar  solution  in  the  equation  and  subtract  the  resulting 
equation  from  (l«  06)  • By  so  doing,  ve  get  an  hoizogeneous  equation* 

The  Particular  Solution 


We  proceed  with  the  derivation  of  the  particular  solution  and  shall 
consider  later  the  solutions  of  the  homogeneous  equation. 

Following  Kelley  and  McShane,  let  us  write  (1*06)  in  the  form 

(1.07)  tj"  + aq'  + br,  ■ c . 

Assume  a solution  of  the  form 


1 


c . 

,r 


on  the  assumption  that  the  q”  and  q*  terms  are  relatively  small*  On  Inserting 
the  solution  in  (1*07)  and  neglecting  r»,f  and  e*.  one  finds 


and 


a / c 

“ b ' b J * 


c a / c v 
b - b (¥  } 


t 


(9) 


In  terms  of  the  notation  of  equation  (l*06)  this  gives 

r *2  , « ti  f r 2 nV 

iAN  n0  * . i(0  +h0*)|  IAN  J iAH  n«  . i(0«  + fcS»)  * 

— [»*  -~+  5 — ^“T [_»*  "T*-1-! 


if  eH  but  the  leading  term  in’a1  and  *b*  of  equation  (1.0?)  are  ignored 
♦Value c of  terms  in  ’a1  and  *b*  for  the  57fflHTH.fi.  Shall  M5^-  at  a velocity 
2000  ft/sec  and  a twist  of  rifling  of  1 in  30: 


/V 

u 

588*5  rad/seo 

K 

1*39  bw'* 

h 

5*51  sec"1 

X 

0*52  sec”3* 

H 

1017  lb. ft2/® 

y* 

52,000  sec-2 

a k 

819  sec  2 

7 

-2.06  secT1, 

a 7 

-1212  sec"2 

21 

O ft  I 

If  n©  , c : and  ht)  in  the  first  bracket  or»  the  right  are  neglected 
aud  the  second  term,  i*  neglected  entirely.  (9)  becomes 


(9.1) 


iAI.  g cos  B 
\i  v 


In  view  of  the  fact  that  the  z aide  is  the  axis  or  lmaglnarles.  (9.1) 
indicates,  on  this  approximation,  that  (a)  the  projectile  is  pointing 
directly  to  the  right,  and  (b)  the  magnitude  of  the  yav  is  given  by 


AH  g cos  8 

(JL  v 


The  preceding  j.s  therefore  a derivation  of  the  result  given  in  Hayes. 
Chapter  X,  p.  teu.  Essentially,  this  expression  for  what  is  called  the 
"yav  of  repose"  vas  used  by  Fovler,  Gallop,  Lock  and  Richmond  in  predicting 
the  "right  drift”  of  a projectile  with  a satisfactory  E^rosaieut  with  experi- 
ment. 


Solutions  for  the  Homogeneous  Equation 

As  vas  pointed  out,  if  we  subtract  the  special  solution  (9)  from  (1.06), 
vs  get  the  following  homogeneous  equation, 

(10)  rj"  - i(n  + is  + ih  + 7v)ti * - |^/B  + ifl{ic  - 7 - 17v)  - hit  + ihx"]r,  = 0 . 

As  in  Chapter  I,  ve  write  (10)  in  the  form 

is"  - iEn * - Ftj  * 0 
and  make  the  transformation 


y exp 


I*  f 


E dt 

< -I 


obtaining 

y”  + ( f - F + | 1 E*)y  - 0 . 

Upon  inserting  the  values  of  E and  F from  (10),  this  becomes 

2 = 

% — B ■*■  § (”iK  + ih  + 2i7  - X)  + | h*  + 

| i(*X  + hX)  - | {x2  + h2  - 7?)  + 

| i(n*  + :1k • + ih'  + X') 


y"  + 


0 . 


If  ve  neglect  the  products  and  squares  of  the  factors  and  also 

their  derivatives  and  write 

a 


ft'  * — 


m 


the  coefficient  of  y in  (ll)  becomes 


s 


T B T 2 
In  view  of  the  fact  that 

u /H  \2 

T ’ 3 “ (2  0)  > 

this  may  be  written 

(12)  (§  of  Q 


X-  - 6 + ^ <H  - .t  + 2r,  - t)  - I X 


1 + nc? 2 (h  - k + 27  - 0 - 


gA 

Oa2 


■ 

J * 


This  quantity  (12)  is  -q^  of  equation  (4.6),  page  6.  Upon  taking  the 
square  root*  of  the  bracket  by  the  aid  of  the  binomial  theorem,  i » e . , 


C— 3 


i 

5 


- 1 + g x + 


and  throwing  away  terms  of  higher  degree  than  the  first  in  x,  we  obtain 

i ( 1 

+ q = + | i -|fio  + i(h-K  + 27-  ()/a  - V<*  J 

i ^£lo  + i(h-K  + 27  - jO/cr  - V°  ^ • 


- o.  = 


In  taking  the  derivative  a . we  consider  only  the  time  derivatives  of 
fl  and  a and  neglect  the  others.  Thus,  as  before  (see  page  7), 


gi  „ i (f  - . 

2q  2 O' 


^Following  a suggestion  of  O.H«  Murphy  (BKL  Technical  Mote  703) , we  may 
replace  the  bracket  of  (12)  by 

(F  ♦ iQ)  » 

If  now  we  place 


a + 


ib  » P + iQ  , 


a2  - b2  ■ P,  2ab  » Q • 

According  to  the  approximation  by  the  binomial  theorem, 

IQ 


while,  according  to  the  above  exact  relations,  the  first  term  should  not  be 


" 

T + t?  • 


but  approximately  "\j~P~  “^1  + r ?2  . On  the  other  hand,  the  second  of 
these  relations  is  fulfilled  by  the  approximation  of  the  binomial  theorem. 
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Kenc*-,  the  solution  for  y is 


(13)  y 


■ ^ Hi 1 f{ 


a<j  - £ + i j (h  - k + 2y  - f)/a  - ( 


r~ 


y_  exD I 

" I 2 

u 


1 ifL  - i[lh  - , t 27  - O/a  + jf  - | j}  at"1 


j(  ' “ 'L 

where  y1  and  y?  are  complex  constants. 


If  ve  write  this  as 


■x 


y a y±e  + y2s  • 


then,  in  view  of  the  fact  that 

tj  “ y exp 


and  that 
we  have  for  r.- 


| i / E dt  | 


° J 

E . (fl  + 1m  ih  + >.), 


UM 


t]  — exp 


2 i J°  (fl  + ih  + iit  + X)  dt  + | + 


Ki,  cap  ^ i J (Q  + ih  + Ik  + X)  dt  + Qg 
o 


u 


where  and  are  complex  constants. 
Finally  we  write 


(1^.1) 


*1  = K^e  + KgO 


From  (1*0  and  (lU.l),  it  is  seen  that 
(15)  Px  = 5 / |l  Fn(l  + o)  + X - - |^h  + k - j(  + 

p,  » I fN  1 f“d  - «>  + X + Vo”!-  fb  + K - ( + c/<7  - (h  - K + 2y  - 


+ (h  - It  + 27  -jf)/c 


s 

i 

>dt, 


J) 


Since  the  stability  factor  is  assumed  to  be  appreciably  greater  than  1, 
the  terms  in  X are  negligible  compared  to  those  in  ft*.  Hereafter,  X will  be 
*E.gt,  for  the  Practice  Projectile  Tll*+  fired  at  a muzzle  velocity  of~ 

3026  ft/sec  with  a twist  of  rifling  of  one  tv.ru  in  25»**  calibers, 

ft  « 1525  rad/sec,. 

X = 0.-2  sec  ^ 

n)* 


( approximately ) . 


neglected.  Equations  (15)  then  become 


1 2 


(15.1) 


1 {[«&♦.).[; 


+ K - / + 6/0  + (h  - K + 


, - 0/3  j 


T>  _ i P i *r\lt  ~\  r x 

*2  ” 2 j ) " u/ 


_ ) 


/ ^ / /„  » i ! 1 

jju+iv-A  t «/g  - - I!  T <7  - X)/® 


UJ 


4t, 


dt 


Tbo  Gene-x  eu.  Solution 


If  the  particular  solution  (9)  is  represented  by  c •>  id  where  c and  d 
are  both  real,  then  the  general  solution  (the  of  the  particular  solution 
and  the  solutions  or  the  homogeneous  equation)  Is 

P F 

1 2 

r(  s e + e + c + id  . 

If  the  initial  value  of  n is  n . etc* : the  constants  K,  end  KU  must 

•o'  ' i 

satisfy  the  conditions: 

i0  * h * h * c0  * i4= 


< - *1  Pi,0  + *2  P2,0  + ci  - K ■ 


Stability  conditions  for  spinning  shell 

Aa  we  have  seen  in  Chapter  I,  the  imaginary  parts  of  and  P2  determine 

the  frequency  or  angular  velocity  of  the  motion  while  the  real  pasts  control 
the  amplitude. 

A shell  Is  ordinarily  said  to  be  dynamically  stable  if  the  amplitudes  of 
both  motions  diminish  with  time.  Thus,  obviously,  the  conditions  for  this 
sort  of  dynamic  stability  are  that  the  real  parts  of  P.^  and  Pg  must  be  nega- 
tive. However,  in  such  a treatment  of  the  dynamical  stability  of  shell,  no 
attention  is  paid  to  the  variation  of  these  factors  with  velocity  ana,  as  a 
•.rule,  the  factors  vary  considerably  as  the  velocity  changes.  It  will  be 
noticed  that  the  amplitude  depends  not  on  the  integrands  involved  in  and 

?2  but  on  the  integrals,  iherefore,  the  integrand  may  be  positive  over 

certain  parte  of  the  velocity  interval  considered  while  the  integral  may  be 
negative.  It  is  apparent  that  the  integrand  itself  should  ba  negagive  at 
the  muzzle  velocity  at  which  the  shell  is  fired  because,  if  it  is  not,  the 
large  amplitudes  might  develop  before  the  integral  becomes  negative  once 
mere.  Thus,  practically  considered,  thejpe  are  two  factors  that  have  to  be 
considered  in  designing  a shell:  (l)  the  integrand  of  the  real  parts  An 

and  P2  should  be  negative  in  the  neighborhood  of  the  muzzle  velocity,  but, (2) 

or, 

■ — ✓ 


on  tbe  other  hand,  if  the  integrands  are  reasonably  negative  in  the  neighbor- 
hood of  the  muzzle  velocity,  it  is  not  essential  that  they  be  negative  at  el  1 
other  velocities  which  the  shell  assumes  on  its  trajectory*  It  is  known 
that  there  are  q'lite  a number  of  shell  which,  in  the  narrow  sense  of  the 
word,  are  dynamically  unstable  on  certain  parts  of  thsi?  trajectories  but 
nevertheiese  are  exceedingly  accurate-  It  would  be  throwing  away  many  perfectly 
good  shell  designs  if  one  insisted  that  the  Integrand  be  negative  s+  all  ve- 
locities which  the  shell  takes  on  during  its  trajectory. 

If  c is  real,  the  four  conditions  for  the  satisfactory  damping  of  the 
yaw  of  a shell  on  a given  trajectory  are: 


near  the  muzzle 


(a)  |^h  + k - i + b/a  + (h  - k + 2y  - ()/a  [‘^5>  0 


(b)  ! h + k - / + a/a  - (h  - k + 2y  - jf)/a |' 


(dynamically  stable) 


along  the  trajectory 


(c)  f |~h  + R 

- / + 

b/a  + (h  « k + 27  - J 0/c  jatlh^O 

(not  necessarily 

-4 

dynamically  stable 

b/a  - (h  - k + 2y  - / ) /o  } dt  ^>0 

at  all  points  on 

(d)  / h + K 

- ( * 

the  trajectory) 

> 

O 1 — 

1 

when  t is  any  time  along  the  trajectory  . 

As  Murphy  suggests  in  BRL  Report  855#  for  practical  purposes,  the  0*6 
on  the  right  of  the  above  inequalities  might  well  be  replaced  by  positive 
constants * 


As  Sterrve  and  Poor  have  pointed  out,  it  may  be  dangerous  to  design  a 
weapon  ior  which  all  four  stability  conditions  are  satisfied  only  vhen 
the  shell  lt>  fired  for  a certain  muzzle  velocity.  If  the  muzzle 
velocity  were  inadvertently  changed,  the  c ondit  ions  'a*  and  *b’  might 
no  longer  be  satisfied. 


Chapter  V 


The  Motion  of  the  Coordinate  System. 
Comparison  with  Fowler,  Gallop,  Lock  and  Richmond. 
Comparison  with  Kelley  and  McShane. 
Neglected  Forces. 
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The  Motion  of  the  Coordinate  System 

As  c 025; area  with  the  coordinate  system  of  Fowler,  Gallop,  Lock  and 
Richmond,  which  is  attached  to  the  particle  trajectory,  oars  has  the 
advantage  that  it  deals  with  only  one  vector,  A.  instead  of  two,  A and  X 
(see  below).  As  compared  with  that  of  Kelley  and  McShane*,  It  has  the 
advantage  of  resembling  more  closely  a Galilean  system  than  one  whose 
coordinates  oscillate  with  the  projectile.  However,  our  coordinate 
system  has  the  disadvantage  that  the  z axis,  originally  assessed  to  be 
horizontal,  will  not,  in  general,  remain  horizontal,  and  the  y axis, 
originally  in  a vertical  plane  through  the  trajectory,  will  not,  in 
general,  remain  in  that  plane.  However,  it  will  be  shown  tnat  the  departure 
of  the  z axis  from  the  horizontal  will  bo  small  and  the  departure  of  the  y 
axis  from  the  vertical  plane  will  be  also. 

Consider  a projectile  moving  vertically  upward.  Let  us  take  the  y axis 
initially  as  coming  out  of  the  plane 

* 


Figure  9 

of  Figure  9 end  the  z axis  as  pointing  to  the  right.  Suppose  there  is  a 
spin  a>^  of  the  coordinate  system  about  the  y axis,  then  the  z axis  will 

turn  toward  the  vertical  with  the  came  angular  velocity  a^..  Place  a vertical 

plans  80  that  it  will  include  the  z axis  end  let  ♦ be  the  angle,  measured  in 
that  plane  between  & horizontal  linn  and  the  z axis.  In  this  case,  $>  = 0^. 

Let  hr  coB?id? r the  general  cases  Suppose  the  tangent  to  the  trajectory 
makes  an  angle  d with  the  horizontal,  let  the  z axis  point  initially  directly 
to  the  right  and  let  the  y axis  be  initially  in  the  vertical  plane  thru  the 
x axis.  See  Figure  10. 


Kelley  and  McShane  have  adopted  the  coordinate  system  of  Nielsen  and  Synge. 
Cf.  "On  the  Motion  of  a Spinning  Shell*  Quarterly  of  Applied  Mathematics, 
Yol.  IV,  No.  3 , October  19^6,  pp»  201-226. 


X 


Kignrn  _!_0 

It  appears  that  the  component  of  a about  an  horizontal  normal  to  the 

y 

initial  direction  of  the  z axis  will  he  ein  0=  Sc,  in  general, 

▼ = sv  sin  f)  . 

It  is  apparent  that  a rotation  © about  the  z axis,  if  z is  horizontal, 
will  not  move  the  y axis  f’.’om  the  vortical. 

On  the  other  hand,  if  the  % axis  is  not  horizontal,  a rotation  about  it, 

® , will  cause  the  y axis  to  leave  the  vertical  plane  through  the  trajectory. 

However,  a is  approximately  proportional  to  © as  long  as  a is  small.  Urns, 

the  departure  of  the  y axis  from  the  vertical  will  be  approximately  proportional 
to  the  product  ®y’®z  sad  will  be  neglected. 

It  is  obvious  that  a rotation  about  the  y axis  will  not  change  the 
direction  of  the  y axis. 

It  thus  appears  that  the  rate  of  departure  of  the  z axis  from  the  hori- 
zontal is  given  by  I » »v  sin  0 , and  that  the  rate  of  dfip  SmTwUTw  Cf  wriu  y axis 

will  he-  cf  higher  order  in  the  spins  of  the  trajectory  ^ 

The  spin  of  the  coordinate  system  caused  by  gravity,  ©^,  has  no  component 
in  a>y.  The  main  component  of  ©y.  is  that  caused  by  the  cross  wind  force,  in 
magnitude,  - kt). 


With  an  assumed  value  of  k u 1.0,  Hitchcock  and  Harrington  have  computed 

. . , v w U dt  Qr  dt 
■ « ” stn  6 ~ cin  6 

..  .■  :-  h 

for  two  cases,  one  where  the  motion  is  undamped  (with  damping  due  to  k neg- 
lected) and  one  where  the  motion  is  damped  with  h « 3*2  and  k w 1.0.  This 
was  for  a 3.3"  shell  with  a velocity  of  2080  ft/sec.  The  results  were  that, 
at  the  cad  of  .6  sec.,  the  maximum  •value  of  » for  the  undamped  case 

sin  6 

was  less  thau  .003  rad,  and  less  than  .002  rad  for  the  damped  case. 
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Fowler  and  associates  refer  the  motion  to  a coordinate  system  attached 
to  the  particle  trajectory  with  only  drag  and  gravity  acting  on  the  particle. 
In  this  system,  they  consider  two  vectors.  One  which  we  Bhail  call  Ap  1b 


the  unit  vector  having  the  direction  of  the  axis  of  the  shell;  the  other,  x. 


a unit  vector  having  the  direction  of 
the  coordinate  system,  mentioned. 


the  actual  trajectory,  both  referred  to 


Fowler  writes 


^F  * ^F*F  + T^F  + aVcF 

a » xip  + yjp  + ahp  • {x,y, z are  direction  cosines.) 

If  we  define  the  vector  yew  indicated  by  ^ as  the  component  of  the 
vector  Ap  normal  to  the  actual  trajectory,  then  with  the  coordinate  system 

of  FGLR, 

« Ap  - X - (»p  - y)JF  + (np  - z)kp 

since  | sag  x 1 . 


On  the  other  hand,  the  vector  representation  of  our  complex  yaw  is 
(raj  + nk) . Hence,  if  t)ie  small  angles  between  our  and  their  coordinate 
axes  may  be  disregarded,  i.e. , if  = J,  kp  «*  k; ;s  = st-  » y,  n «*  np  - z . 

FGLR  introduce  two  complex  variables  in  their  treatment,  q and  £ . The 
definitions  of  these  variables  are  as  follows; 

Tj+cg-Mp+inj,  , c 5 * y + iz 

where  c = cos  Q. 

From  these,  it  is  seen  that 

tj  - (nip  - y)  + ifiip  - z)  , 

and  Fowler  ■ s q is  equivalent  to  our  q,  except  for  the  slight  difference 
between  our  axes  and  his. 


Cn  page  337,  FGLR  give,  with  slight  changes  in  notation,  the  following 
equations  involving  q and  £ . 

/ 2 ) 

+ifi(  K 'lA-  y) ’ + IV  !- 


(3-613)  q"  -(ifl-h-ie+iXjq' 

- 1 ific 1 - he*  - cr j £ 53  i*»v 

* Where  0^  is  the  Tnclination  of  the  particle- 'trejectoryv 


; '■ 


*• 

m ■ 

E- 


Tbis  is  to  be  compered  with  car  equation  (1.06)  on  page  20.  The  difference?5 
are  as  follows; 

a.  FGLR  have  two  terms  in  r,;  ~{k'  - )*n  end  + (k  - i?y)  . 

c 

b.  They  have  a £ term  which  they  later  discard  as  being  of  the 


c.  On  the  right  hand  side,  ve  have  two  terns  in  the  bracket  of  (1.06) 


n&’~ 

C 


ana 


ih0' 

"IT" 


which  Fowler  considers  neRlifdble, 


FGLR  give  as  the  solution  (p.  339)  for  the  simplified  homogeneous 
equation  in  n.  with  a alight  change  of  form, 

S'*  vVar  r0-| 


/ n 1.  ^ 

» n « tjsr  » 


/'fig  Y 

im:J 

i n o # 


+ 


P . p 
1'  “2 


4.  / ^ 
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and  1C  and  Kg  are  arbitrary  constants. 
It  may  be  shown  that 
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In  view  of  this,  (3*6234)  may  be  written 

P p 

1 ~ 2 

n « K Le  A + Kge  * 

with  the  P's  now  defined  by 

Pl'  P2  * § ♦ *)  -|^+K.jf  + |+(h-K  + 27  - jf)/<T]  t dt  . 

This  is  to  be  compared  with  (15.I)  on  page  25.  It  will  be  seen  that 
the  results  are  identical. 
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Another  comment  on  Fowler's  treatment;  he  claims  his  result  to  be 
valid  only  for  large  values  of  fij  our  results  indicate  that  his  results 
are  of  greater  generality  than  he  realized.  In  fact,  they  merely  depend 
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on  the  usual  W. K.B.  approximations  (see  p.  q)  <. 


Comparison  with  Kelley-McShane 

Kelley-McShane  use  the  arc  length  as  the  Independent  variable,  instead 
of  the  time*  Dr.  Galbraith  has  kindly  transferred  ifclley-McShsne 1 8 sq-nanion 
(1.25),  BRL  Report  4h6,  to  its  equivalent  with  time  aa  the  independent  variable 
The  result,  with  our  notation  for  the  damping  factors,  is: 

Tj  - iq(n  + itt  + ih  + Tv  + * ) - + - 7 - IX)  - i -flfe [ 

V r 

- ia  0'  (l  + | + ^)  - 0”  . 


*ln  view  of  the  fact  that  FGLR  made  their  observations  by  yaw  cards,  they  give, 
in  addition  to  an  expression  for  tj,tho  complex  yaw,  expressions  for  6^  and  for 
<t>,  the  angle  between  a vertical  plane  through  the  trajectory  and  the  plane  of 

the  yaw:  2 2 2 2 2 

5 ’ = or  sin  qt  + P cos“  qt  , 

, r i 

♦ = ♦ + % Dt  + arc  tan  J coth  (.1  - a-)  tan  ot  ► . 
u * V my  ' i 


In  these. 


a **  maximum  yaw 
0 = minimum  yaw 
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If  'cLis  is  compand  with  (l*C *6)  os  page  20,  it  is  apparent  that,  aside  from 
the  terms  in  t which  we  neglect,  the  only  difference  is  the  term 

-"S  sJja  9 in  the  coefficient  of  n.  In  view  of  the  fact  that  this  term  is 
u 1 

omitted  in  ’'Exterior  Ballistics”,  by  McEtap? ,1 Kelley  and  Reno,  it  is  apparent 
that  Kelley  end  McShnne  consider  it  ohcuxn  ue  Cuu^^oq* 

Keglscted  Forces 


We  have  so  far  neglected  three  elements  of  the  fores  system  because, 
in  practice,  the  corresponding  damping  factors  are  negligible*  These  are: 


the  Magnus  force  due  to  yawing, 


the  Magnus  moment  due  to  yawing. 


a.  The  force  due  to  yawing 

The  yawing  moment  due  to  yawing  Is  produced  by  the  product  of 
a force  multiplied  into  the  distance  between  the  center  of  application  of 
the  force  end  the  center  of  gravity.  This  force  is  the  force  due  to  yawing. 
It  is  Important  in  the  motion  of  dirigibles  and  of  interest  in  developing  a 
tneory  of  the  yawing  moment  due  to  yawing.  It  has  a negligible  effect  on 
the  motion  of  ordinary  shell. 


b.  The  Magnus  fores  due  to  yawing 

Just  as  there  is  a Magnus  force  which  corresponds  to  the  cross 
wind  force,  so  there  is  also  a Magnus  force  due  to  yawing.  As  a result  of 
the  spin  of  a shell,  the  Magnus  force  due  to  yawing  differs  in  magnitude  and 
direction  frum  the  force  due  to  yawing.  This  force  is  also  negligibly  small 
and  need  not  be  considered  in  the  theory  of  daaqping. 

c.  The  Magnus  moment  due  to  yawing 

Corresponding  to  the  yawing  moment  due  to  yawing,  if  the  shell 
is  spun,  there  will  be  a Magnus  moment  due  to  yawing  but  this  Magnus  moment 
is  also  negligibly  small. 
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Miss  M.  E.  Harrington  also  c -imputed  the  motion  of  my  coordinate  system 
to  determine  the  departure  from  the  horizontal  of  an  axis  originally 
horizontal* 
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